Computational method to reduce the search space for directed protein evolution

We introduce a computational method to optimize the in vitro evolution of proteins. Simulating evolution with a simple model that statistically describes the fitness landscape, we find that beneficial mutations tend to occur at amino acid positions that are tolerant to substitutions, in the limit of small libraries and low mutation rates. We transform this observation into a design strategy by applying mean-field theory to a structure-based computational model to calculate each residue's structural tolerance. Thermostabilizing and activity-increasing mutations accumulated during the experimental directed evolution of subtilisin E and T4 lysozyme are strongly directed to sites identified by using this computational approach. This method can be used to predict positions where mutations are likely to lead to improvement of specific protein properties.

Prion protein: Evolution caught en route

The prion protein displays a unique structural ambiguity in that it can adopt multiple stable conformations under physiological conditions. In our view, this puzzling feature resulted from a sudden environmental change in evolution when the prion, previously an integral membrane protein, got expelled into the extracellular space. Analysis of known vertebrate prions unveils a primordial transmembrane protein encrypted in their sequence, underlying this relocalization hypothesis. Apparently, the time elapsed since this event was insufficient to create a “minimally frustrated” sequence in the new milieu, probably due to the functional constraints set by the importance of the very flexibility that was created in the relocalization. This scenario may explain why, in a structural sense, the prion protein is still en route toward becoming a foldable globular protein.

Protein evolution on partially correlated landscapes.

We extend an earlier model of protein evolution on a rugged landscape to the case in which the landscape exhibits a variable degree of correlation (i.e., smoothness). Correlation is introduced by assuming that a protein is composed of a set of independent blocks or domains and that mutation in one block affects the contribution of that block alone to the overall fitness of the protein. We study the statistical structure of such landscapes and apply our theory to the evolution by somatic hypermutation of antibody molecules composed of framework and complementarity-determining regions. We predict the expected number of replacement mutations in each region.

Positive Darwinian selection drives the evolution of several female reproductive proteins in mammals

Rapid evolution driven by positive Darwinian selection is a recurrent theme in male reproductive protein evolution. In contrast, positive selection has never been demonstrated for female reproductive proteins. Here, we perform phylogeny-based tests on three female mammalian fertilization proteins and demonstrate positive selection promoting their divergence. Two of these female fertilization proteins, the zona pellucida glycoproteins ZP2 and ZP3, are part of the mammalian egg coat. Several sites identified in ZP3 as likely to be under positive selection are located in a region previously demonstrated to be involved in species-specific sperm-egg interaction, suggesting the selective pressure is related to male-female interaction. The results provide long-sought evidence for two evolutionary hypotheses: sperm competition and sexual conflict.

The crystal structure of a heptameric archaeal Sm protein: Implications for the eukaryotic snRNP core

Sm proteins form the core of small nuclear ribonucleoprotein particles (snRNPs), making them key components of several mRNA-processing assemblies, including the spliceosome. We report the 1.75-Å crystal structure of SmAP, an Sm-like archaeal protein that forms a heptameric ring perforated by a cationic pore. In addition to providing direct evidence for such an assembly in eukaryotic snRNPs, this structure (i) shows that SmAP homodimers are structurally similar to human Sm heterodimers, (ii) supports a gene duplication model of Sm protein evolution, and (iii) offers a model of SmAP bound to single-stranded RNA (ssRNA) that explains Sm binding-site specificity. The pronounced electrostatic asymmetry of the SmAP surface imparts directionality to putative SmAP–RNA interactions.

A natural polymorphism in β-lactamase is a global suppressor

A M182T substitution was discovered as a second-site suppressor of a missense mutation in TEM-1 β-lactamase. The combination of the M182T substitution with other substitutions in the enzyme indicates the M182T substitution is a global suppressor of missense mutations in β-lactamase. The M182T substitution also is found in natural variants of TEM-1 β-lactamase with altered substrate specificity that have evolved in response to antibiotic therapy. The M182T substitution may have been selected in natural isolates as a suppressor of folding or stability defects resulting from mutations associated with drug resistance. This pathway of protein evolution may occur in other targets of antimicrobial drugs such as the HIV protease.

Why are some proteins structures so common?

Many biological proteins are observed to fold into one of a limited number of structural motifs. By considering the requirements imposed on proteins by their need to fold rapidly, and the ease with which such requirements can be fulfilled as a function of the native structure, we can explain why certain structures are repeatedly observed among proteins with negligible sequence similarity. This work has implications for the understanding of protein sequence structure relationships as well as protein evolution.

Genetic definition of a protein-splicing domain: Functional mini-inteins support structure predictions and a model for intein evolution

Inteins are protein-splicing elements, most of which contain conserved sequence blocks that define a family of homing endonucleases. Like group I introns that encode such endonucleases, inteins are mobile genetic elements. Recent crystallography and computer modeling studies suggest that inteins consist of two structural domains that correspond to the endonuclease and the protein-splicing elements. To determine whether the bipartite structure of inteins is mirrored by the functional independence of the protein-splicing domain, the entire endonuclease component was deleted from the Mycobacterium tuberculosis recA intein. Guided by computer modeling studies, and taking advantage of genetic systems designed to monitor intein function, the 440-aa Mtu recA intein was reduced to a functional mini-intein of 137 aa. The accuracy of splicing of several mini-inteins was verified. This work not only substantiates structure predictions for intein function but also supports the hypothesis that, like group I introns, mobile inteins arose by an endonuclease gene invading a sequence encoding a small, functional splicing element.

On roads not taken in the evolution of protein catalysts: Antibody steroid isomerases that use an enamine mechanism

Reactive immunization has emerged as a new tool for the study of biological catalysis. A powerful application resulted in catalytic antibodies that use an enamine mechanism akin to that used by the class I aldolases. With regard to the evolution of enzyme mechanisms, we investigated the utility of an enamine pathway for the allylic rearrangement exemplified by Δ5-3-ketosteroid isomerase (KSI; EC 5.3.3.1). Our aldolase antibodies were found to catalyze the isomerization of both steroid model compounds and steroids. The kinetic and chemical studies showed that the antibodies afforded rate accelerations up to a factor of 104 by means of an enamine mechanism in which imine formation was the rate-determining step. In light of our observations and the enzyme studies by other workers, we suggest that an enamine pathway could have been an early, viable KSI mechanism. Although this pathway is amenable to optimization for increased catalytic power, it appears that certain factors precluded its evolution in known KSI enzymes.

Chloroplast small heat shock proteins: Evidence for atypical evolution of an organelle-localized protein

Knowledge of the origin and evolution of gene families is critical to our understanding of the evolution of protein function. To gain a detailed understanding of the evolution of the small heat shock proteins (sHSPs) in plants, we have examined the evolutionary history of the chloroplast (CP)-localized sHSPs. Previously, these nuclear-encoded CP proteins had been identified only from angiosperms. This study reveals the presence of the CP sHSPs in a moss, Funaria hygrometrica. Two clones for CP sHSPs were isolated from a F. hygrometrica heat shock cDNA library that represent two distinct CP sHSP genes. Our analysis of the CP sHSPs reveals unexpected evolutionary relationships and patterns of sequence conservation. Phylogenetic analysis of the CP sHSPs with other plant CP sHSPs and eukaryotic, archaeal, and bacterial sHSPs shows that the CP sHSPs are not closely related to the cyanobacterial sHSPs. Thus, they most likely evolved via gene duplication from a nuclear-encoded cytosolic sHSP and not via gene transfer from the CP endosymbiont. Previous sequence analysis had shown that all angiosperm CP sHSPs possess a methionine-rich region in the N-terminal domain. The primary sequence of this region is not highly conserved in the F. hygrometrica CP sHSPs. This lack of sequence conservation indicates that sometime in land plant evolution, after the divergence of mosses from the common ancestor of angiosperms but before the monocot–dicot divergence, there was a change in the selective constraints acting on the CP sHSPs.

Identification of interdomain sequences promoting the intronless evolution of a bacterial protein family.

In the evolution of eukaryotic genes, introns are believed to have played a major role in increasing the probability of favorable duplication events, chance recombinations, and exon shuffling resulting in functional hybrid proteins. As a rule, prokaryotic genes lack introns, and the examples of prokaryotic introns described do not seem to have contributed to gene evolution by exon shuffling. Still, certain protein families in modern bacteria evolve rapidly by recombination of genes, duplication of functional domains, and as shown for protein PAB of the anaerobic bacterial species Peptostreptococcus magnus, by the shuffling of an albumin-binding protein module from group C and G streptococci. Characterization of a protein PAB-related gene in a P. magnus strain with less albumin-binding activity revealed that the shuffled module was missing. Based on this fact and observations made when comparing gene sequences of this family of bacterial surface proteins interacting with albumin and/or immunoglobulin, a model is presented that can explain how this rapid intronless evolution takes place. A new kind of genetic element is introduced: the recer sequence promoting interdomain, in frame recombination and acting as a structure-less flexibility-promoting spacer in the corresponding protein. The data presented also suggest that antibiotics could represent the selective pressure behind the shuffling of protein modules in P. magnus, a member of the indigenous bacterial flora.

A protein-binding domain, EH, identified in the receptor tyrosine kinase substrate Eps15 and conserved in evolution.

In this report we structurally and functionally define a binding domain that is involved in protein association and that we have designated EH (for Eps15 homology domain). This domain was identified in the tyrosine kinase substrate Eps15 on the basis of regional conservation with several heterogeneous proteins of yeast and nematode. The EH domain spans about 70 amino acids and shows approximately 60% overall amino acid conservation. We demonstrated the ability of the EH domain to specifically bind cytosolic proteins in normal and malignant cells of mesenchymal, epithelial, and hematopoietic origin. These observations prompted our search for additional EH-containing proteins in mammalian cells. Using an EH domain-specific probe derived from the eps15 cDNA, we cloned and characterized a cDNA encoding an EH-containing protein with overall similarity to Eps15; we designated this protein Eps15r (for Eps15-related). Structural comparison of Eps15 and Eps15r defines a family of signal transducers possessing extensive networking abilities including EH-mediated binding and association with Src homology 3-containing proteins.

Accelerated evolution in the protein-coding regions is universal in crotalinae snake venom gland phospholipase A2 isozyme genes.

The nucleotide sequences of four genes encoding Trimeresurus gramineus (green habu snake, crotalinae) venom gland phospholipase A2 (PLA2; phosphatidylcholine 2-acylhydrolase, EC 3.1.1.4) isozymes were compared internally and externally with those of six genes encoding Trimeresurus flavoviridis (habu snake, crotalinae) venom gland PLA2 isozymes. The numbers of nucleotide substitutions per site (KN) for the noncoding regions including introns were one-third to one-eighth of the numbers of nucleotide substitutions per synonymous site (KS) for the protein-coding regions of exons, indicating that the noncoding regions are much more conserved than the protein-coding regions. The KN values for the introns were found to be nearly equivalent to those of introns of T. gramineus and T. flavoviridis TATA box-binding protein genes, which are assumed to be a general (nonvenomous) gene. Thus, it is evident that the introns of venom gland PLA2 isozyme genes have evolved at a similar rate to those of nonvenomous genes. The numbers of nucleotide substitutions per nonsynonymous site (KA) were close to or larger than the KS values for the protein-coding regions in venom gland PLA2 isozyme genes. All of the data combined reveal that Darwinian-type accelerated evolution has universally occurred only in the protein-coding regions of crotalinae snake venom PLA2 isozyme genes.